


Institutional Archive of the Naval Postgraduate School 


Calhoun: The NPS Institutional Archive 
DSpace Repository 


Theses and Dissertations 1. Thesis and Dissertation Collection, all items 


1964 


Thermal convection in Lake Washington, 
winter 1962-1963. 


Syck, James Marvin. 


University of Washington 


http://ndl.handle.net/10945/13279 


Downloaded from NPS Archive: Calhoun 


| Calhoun is the Naval Postgraduate School's public access digital repository for 
D U DLEY research materials and institutional publications created by the NPS community. 
get Calhoun is named for Professor of Mathematics Guy K. Calhoun, NPS'‘s first 
KNOX appointed — and published — scholarly author. 





LIBRARY Dudley Knox Library / Naval Postgraduate School 
411 Dyer Road / 1 University Circle 


http://www.nps.edu/library Monterey, California USA 93943 


NPS ARCHIVE 


1964 
SYCK, J. 





THERMAL CONVECTION IN LAKE WASHINGTON, 
WINTER 1962-1963 


JAMES MARVIN SYCK 











i = 








%)) J}) SPOOL) 


THERMAL CONVECTION IN LAKE WASHINGTON, 


WINTER 1962-1963 


by 


JAMES MARVIN SYCK 


A thesis submitted in partial fulfillment 


of the requirements for the degree of 


MASTER OF SCIENCE 


UNIVERSITY OF WASHINGTON 
196), 


Approved by 
Department in which degree is granted 


Date 








)) » 


Library 


U.S. Naval Postaesduate School 


Monterey, California 


ACKNOWLEDGMENTS 


The author wishes to thank the members of the faculty of the 
Department of Oceanography for their assistance in this work. Parti- 
cular thanks to Dr. Clifford A. Barnes for pointing out the possibili- 
ties of using Lake Washington for studying convective flows. Thanks 
also are given to Drs. Maurice Rattray, Jr. and YT. Saunders i“nglish 
for their many helpful suggestions in the preparation of the manuscript. 

Incoming solar radiation data was obtained from the Department 
of Atmospheric Sciences at the University of Washington. 

This research was performed while the author was on active duty 
with the United States Navy in the Junior Line Officers Advanced 
Scientific Program administered by the Bureau of Naval Personnel. 
Financial support of the field work came in part from Navy contracts 
peer n7 710) and Nonr-4,77( 37). 





INTRODUCTION . 


Purpose 


e 


»») 


Cd e 


asl 


TABLE OF CONTENTS 


Nacuresot the Proplem, ..39m. 2. <<. 4am 


Previous Work 


Description of Lake Washington... . 


Materials and Methods 


INVESTIGATION . 


The Temperature Structure ..... 


Heat Budget 


Downslope Current off Juanita Bay .. 


MaSCUSSION ... 


CONCLUSIONS . . 


LITERATURE CITED 


e 


@ 


)) 


CO ON OO fF WW 


1O 


one 
- 1h 
ee 





10. 


1 


le. 


13. 


»> 


pel 
LIST OF FIGURES 
The bottom topography of Lake washington and the location 
Br themcectiOnsS . . «as #944 8 4 0 «© 6 ee sw eee 


The vertical temperature structure off Fairweather Point, 


memmeerol 1962-603. . << 5 SW «© * 6 © 6 gee 6 sc 


Temperature profile from Sand Point 
11 and 12 January 1963 ...... 
Temperature profile from Sand Point 
1h and 21 January 1963 ..... 
Temperature profile from Sand Point 
25 and 27 January 1963 

Temperature profile from Sand Point 
29 and 30 January 1963 .....-. 
Temperature profile from Sand Point 
31 January and 1 February 1963 .. 
Temperature profile from Sand Point 
18 February and 6 March 1963 ... 
Temperature profile from Sand Point 


AW Neves O08. rr rn) ae 


Temperature profile from Kenmore to 


to Juanita Bay, 


to Juanita Bay, 


to Juanita Bay, 
JS ty, 
» ca 
oe 


to Juanita Bay, 


Fairweather Point, 


Cimdanuawyel yon. « . oeeis « « « « oes «© « «eareeeee 


Temperature profile along the axis of the lake from 


Sand Point to Leschi Park, 21 January 1963 .... 


Daily average air temperature and daily average wind 


Bpecdeat oand POlnung =. . . . 


Heat loss computed from the heat budget and the observed heat 


and computed current in the Sand Point-Juanita Bay section 


» 29 


380) 
loss 


eek 





II. 


Pil . 


») )) 


iv 


LIST OF TABLES 


Average daily wind for the warming cycle, 2 February to 
Pimms) 1963 . 3 6 se ee ee sw Gun Sed oss (ee 
Meteorological data used in the heat budget calculations .. . 17 


Tecms of the heat budget equations. . .... .. . . see dC 





») )) 


ABSTRACT 


The winter convection regime in Lake Washington was investigated, 
especially the flow along the bottom off the mouths of some of the many 
Shallow bays. The density of fresh water at a given depth, or pressure, 
is a function of its temperature, so temperature was measured in detail 
in a portion of the lake during the winter cooling period of 1962-63. 
The isotherms at the edges of the lake were found to be closely packed 
and roughly parallel to the sloping bottom. This structure frequently 
extended to the middle and occasionally across the lake whereas most of 
the water column in the center of the lake was nearly isothermal. Volume 
transport was computed from consideration of the amount of heat removed 
from the deep water of the lake. From the volume transport off Juanita 
Bay the velocity of the water down the side of the lake was estimated to 


be 2 cm/sec. 





») )) 


INTRODUCTION 


PURPOSE 

The purpose of this study is to describe convection in Lake 
Washington associated with winter cooling under conditions of low wind 
and large cooling rates. It was necessary to (1) define the vertical 
and horizontal temperature structure of the lake, (2) compute the heat 
budget of the lake during the period of winter cooling, (3) compute the 
vertical transport and associated velocities required for bottom water 


replacement by convection currents. 


NATURE OF THE PROBLEM 

The density structure of Lake Washington gives a measure of the 
potential for flow of the system, Since the density of fresh water at 
a given depth is a function of temperature and the amount of dissolved 
or suspended material, measurement of these variables at successive 
times will give the potential for flow of the system. In addition, it 
will be possible to define the advective flow in the layers of the lake 
below the depth of direct surface influence, as all changes of heat at 
depth are the result of advection (diffusion is considered to be 


relatively unimportant because the water is nearly isothermal). 


PREVIOUS WORK 

The term "overturn" was introduced by Whipple in 1895 (Birge, 
190) to describe the annual replacement of bottom water in lakes 
associated with winter cooling. This choice of terminology appears to 
be unfortunate because it implies that individual water columns becom? 


unstable and capsize. Early workers, particularly those working in 
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Shallow lakes in areas where winter cooling is so great that the lakes 
freeze, tended to discount the importance of convection. Birge (1910) 
went so far as to say, 
"convection currents are far less important agents for distri- 
buting heat than are mechanical currents caused by the wind. 
Indeed, it would be difficult to show that convection currents 
have any such efficiency in carrying heat as to make them wortn 
Serious consideration." 
More recently, Mortimer (1955) recognized a problem with the term 
"overturn" and was careful to state, "'overturn' connotes the fact that 
bottom layers previously isolated below the thermocline are brought to 
the surface." Mortimer goes on to say that such an indefinite statement 
is necessary because so little attention has been given to water move- 
ments leading to the destruction of the thermocline. 

The limnology of Lake Washington was discussed by Scheffer and 
Robinson (1939). They found that bottom water replacement takes place 
over an extended period in the winter, from December to February or 
March. The temperature in Lake Washington does not go below the 
temperature of maximum density of fresh water, 4° C, so that there is 
no winter stable period. 

Rattray, Seckel and Barnes (1954) studied the salt balance in 
the Lake washington Ship Canal. They repeatedly occupied a station in 
the deep part of the lake off Madison Park. They found that there was 
bottom water replacement even in the presence of a density gradient. 
They showed that with a density increase of 1.3 x 1074 g/em from top 
to bottom, it was still possible to have an "overturn". Tne energy for 
the mixing presumably een from the wind. They also considered the 
effect of river inflow and concluded that 25 percent of the lake volume 


was flushed each year. 
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Gould and Budinger (1958) proposed that convection currents were 
responsible for the W-shaped bottom profile of Lake washington. They 
postulated that currents arising in the shallow bays around the edges 
of the lake were strong enough to erode the foot of the slopes, or at 
least to prevent deposition of sediments there. This is an alternative 
mechanism for the replacement of bottom water in the presence of a 
density gradient. The water need only become very dense locally in the 
Shallow bays; it can then flow under the water in the central portion of 
the lake, displacing the local water, with its associated density 
structure, upward into regions where the structure can be altered by 
wind action and heat exchange across the surface. This process will 
continue as long as water in the bays is colder than the main body of 
the lake, which agrees with the observation of Scheffer and Robinson 
that bottom water replacement occurs over an extended period. 

Church (1942) worked on the annual temperature cycle in Lake 
ehigan. He found that with a density increase of 1.5 x 10-4 g/cm, 
over the depth range encountered in Lake Washington, stability was 
maintained against a wind of Beaufort force 2. In Lake Michigan, 
however, the wind might have been the dominant factor in the breakdown 
of the thermocline. Winter winds in that area were very strong and 


there were few shallow bays to serve as source areas of cold bottom water. 


DESCRIPTION OF LAKE WASHINGTON 
Lake Washington is a narrow lake with its long axis in a north- 


south direction. It is about 1.5 miles wide and 18.5 miles long. Along 
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mach of the length of the lake there is a rise of 2 to 10 m, giving 
transverse profiles a subdued W-shape (fig. 1). The margin of the lake 
is indented by many shallow bays. The lake communicates with the salt 
water in Puget Sound by an 18-mile canal with a lock at the seaward end. 
A more extensive discussion of the morphology of the lake is given in 


Gould and Budinger (1958). 


MATERTALS AND METHODS 

Temperature was measured with a thermistor which formed one leg 
of a Wheatstone bridge. The bridge output, which was nearly linear with 
temperature, was fed into a millivolt recorder. The correction to the 
thermistor output for heating by the puildee comment was always less 
than 0.01° C3; the correction was considered negligible. 

In normal operation, the full scale recorder range was 2° C, giving 
a sensitivity of 0.01° C. The system was calibrated by the Applied Physics 
Laboratory, University of Washington. The accuracy of the system was 
checked by field comparisons with reversing thermometers. A bias of 
0.139 C was found. All temperature determinations have been corrected 
for this bias which did not change during the program. 

The standard deviation of the difference between temperatures 
determined by the reversing thermometers and the thermistor was 0.031° C. 
Part of the variability is associated with the limits of accuracy of the 
reversing thermometer and part with the inaccuracy of the thermistor. 

The standard deviation of the averaged determinations by a pair of 
reversing thermometers has been given as 0.012° C (Fofonoff, 1963). 


Since the variance of a sum of independent random variables equals the 





sum of the variances, the standard deviation of the thermistor approxi- 
mates 0.0289 C. The precision of the thermistor (as defined by twice 
the standard deviation) is therefore about 0.05° C, which is about one- 
half the precision of the reversing thermometers. 

Chlorinity was determined by the Volhard method. The carbonates 
and sulfides in the lake water interfere with the Mohr titration which 
is usually used on sea water. 

Positions were determined either by horizontal sextant angles or 


by magnetic bearings from a pelorus on the flying bridge. 





INVESTIGATION 


THE TEMPERATURE STRUCTURE 

The vertical temperature structure at a station on the east side 
of Lake Washington off Fairweather Point was monitored during the winter 
of 1962-63 (fig. 2). Cooling continued from the first observation on 
26 December until 31 January. During this time, no unstable water columns 
were observed, with the possible exception of the last few meters on 
29 and 30 January. A temperature increase in the bottom layers existed 
in previous years (Rattray, Seckel and Barnes, 195) when dissolved salt 
stabilized the water column. 

If the water columns observed on 29 and 30 January were stabilized 
by dissolved salt, there would have to be a chlorinity difference of 
0.007 per mille in the bottom 3 m. The temperature inversion was too 
close to the bottom to permit sampling there. The inversion also could 
have been caused by a turbidity current moving out into the deep water. 
The column would remain stable as the sediment settled and the inversion 
disappeared. Measurements with a hydrophotometer in shallow water indi- 
cated that there was suspended material near the bottom. 

The lake off Fairweather Point was isothermal on 26 December and 
10 January. On 1 January the surface layer had cooled. but there was 
& large temperature decrease in the last 7 to 8 m. The decrease was 
present on 21 January, while the thermocline was much shallower than it 
was earlier. Surface warming began on 1 February, but deep water 
continued to cool. Cooling was probably caused by the arrival of cold 


water which was in transit when surface warming began. By 6 March, 
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warming had taken place at all levels and thermal stratification began. 
Stratification was more pronounced by 2) April and the temperature 
had increased throughout the water column, suggesting an appreciable 
movement of water from the surface to depth during the initial warming. 
The process of replacement of bottom water during warming is different 
from replacement during cooling. The source of energy for replacement ° 
during warming is the wind. The spring storms cause some mixing to the 
bottom. In 1963, winds stronger than any during cooling period occurred 
on 3 and 26 February, 15, 18, 28, 29, 30 March, 15 and 22 April (table I). 
Thirteen sections were made from Sand Point to Juanita Bay (fig. 
3 to 9). They covered the two major cold periods in January, and the 
return of stratification with warming. In all sections during the 
cooling period, there was a large tilt to the isotherms near Juanita 
Bay. The tilted isotherms and the low temperatures in the bay was an 
indication of flow out of the bay along the bottom. 
The cold water often found in the deep on the Sand Point side 
of the section did not appear to have originated in Juanita Bay, as the 
isotherms were not continuous from the mouth of the bay. That water was 
too cold to have come from Juanita Bay. Water from that source would 
have been warmer because of mixing currents. The surface water near 
Sand Point was never as cold as the water in the deep paths of the 
section. The coldest water was in the thalweg, the line of greatest 
depth along the axis of the lake. The most probabie origin of that water 
was in the northern portion of the lake. 
A section from Kenmore to Fairweather Point was prepared to trace 


the origin of the cold water in the thalweg (fig. 10). The section is a 
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composite of sections from Kenmore to Sand Point on 27 January and from 
Sand Point to Fairweather Point on 25 January. The temperatures in the 
portion of the earlier section were adjusted to the levels prevailing 
2? January. 

There was a marked similarity between the temperature structures 
in the north end of the lake and the Sand Point-Juanita Bay sections. 
This was expected because the north end of the lake has extensive shallow 
areas for the formation of very dense water. The cold water at the foot 
of the slope in the Kenmore region presumably was not isolated from the 
cold water in shallower depths, as continuity was maintained along the 
thalweg. 

A section from Sand Point to the deepest part of the lake, off 
Leschi Park, was taken on 21 January (fig. 11). The cold water in the 
central region corresponded in position to both Fairweather Bay on the 
east and to Union Bay on the west and probably represents the cold water 


flowing into the deepest part of the lake from each side. 


HEAT BUDGET 

Knowledge of the heat budget made it possible to compare the 
observed heat loss in a portion of the lake with the computed heat loss 
expected from meteorological conditions. If the observed heat loss and 
the computed heat loss agreed, then all of the heat was removed through 
the surface. If the observed and computed heat loss disagreed, then 
there must have been exchange of water and heat across the lateral 


boundaries of that part of the lake. 
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The terms in the heat budget equation were incoming radiation (Q,), 
reflected radiation (Q;), back radiation (Q,), evaporation (Q_.), and 
conduction of sensible heat (Q,,). The incoming radiation was measured 
with an Eppley pyrheliometer on top of a three-story building (Johnson 
Hall on the University of Washington campus) about 1.5 miles from the 
lake. The remaining terms were calculated from expressions given by 
T. Laevastu (196C), modified slightly for fresh water. 

Qr * 0.15 Qs - (0.01 Q,)% 
Qe = (0.26 + 0.077 V) (e, - eg) (0.2 L) 
Q = 20.65 (14.38 - 0.09 T, - 0.046 r) (1 - 0.0765 C) 
Qn = 36 (0.26 + 0.077 V) (T, - Ty) 
Where all of the Q's are in g cal/cm@/day: 

T,y is the temperature of the water surface in °C, 

T, is the temperature of the air in °C, 

r is the relative humidity in percent, 

C is the cloud cover in tenths, 

V is the wind speed in m/sec, 


e,, is the saturated vapor pressure at the temperature of the 
surface in mb, 


€, is the vapor pressure of the air above the water surface in 
mb, 


L is the latent heat of vaporization in cal/gm. 
Laevastu uses 39, instead of 36, as the coefficient in his equation for 
the transfer of sensible heat. The value comes from the choice of the 
constant in the Bowen ratio. Limnologists have used a slightly different 
value of this constant, resulting in a smaller coefficient in the equation 


of sensible heat transfer (Anderson, 195h). 
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The meteorological data used to make the heat budget calculations 
were obtained from the weather facility at Sand Point Naval Air Station 
(table II). The anemometer was 10 m above the lake level, the wet and 
dry bulb thermometers were 8 m above lake level. The two major cooling 
periods in January 1963, with daily average temperatures less than 0° C, 
had an associated increase in wind speed (fig. 12). The heat budget was 


computed from 10 January to 1 February (table III). 


DOWNSLOPE CURRENT OFF JUANITA BAY 

There are two main assumptions in the calculation of downslope 
current off Juanita Bay. The first is that there is no flow into or 
out of the Sand Point-Juanita Bay section, so that all of the heat loss 
must take place at the surface. The second is that the influence of 
surface processes does not extend below 20 m. If this assumption is 
violated, the computed current will be too high. For the purposes of 
this calculation, the Sand Point to Juanita Bay section was divided into 
three parts. 


Sand Point Juanita Bay | 
— 4 I I . 7 
‘ Tiel a _ 20 m 


41 
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The volume transport was determined by the amount of water at the average 
temperature of Juanita Bay, region II, which would have to mix with the 
water in region I to produce the observed heat loss there. The formula 


used was: 


2 Heat 


." (> Tir) / Cp At 


where V is the volume transport in cm/sec; 





ae 
‘Heat is the change in heat of the water in region I from one 
time to the next, 

Ty is the average temperature of region I, 

Try is the average temperature of region ITI, 

and Cy are density and specific heat (assumed to be 1), 

.t is the time in seconds between sections. 

The current velocity was found by dividing the volume transport 
by the area through which the current flowed. The vertical temperature 
structure in the slope region off Juanita Bay showed that the flow varied 
in thickness from 1 to 5m. Since the thickness of the flow was usually 
closer to 5 m, that value was used in the calculation. This choice gave 


a conservative estimate of the computed current velocity. 





DISCUSSION 


Although dieeeiaed salt contributed to the density structure of 
the lake in the past (Rattray, Seckel and Barnes, 195) 20 chlorinity 
determinations showed that any contribution to the density structure was 
small. The maximum chlorinity increase from the surface to the bottom 
at four deep stations was 0.003. This would cause a change in the density 
equivalent to a temperature change of 0.19 C. One station off Leschi 
Park on 18 February had a maximum chlorinity difference in the water 
column of 0.009 per mille. The chlorinity was never greater than 0.011 
per mille;the average was 0.00) per mille. 

The observed heat loss in the Sand Point-Juanita Bay section was 
compared to the calculated heat budget (fig. 13). The observed heat loss 
in the section was about 20 percent smaller than the computed heat loss. 
There are three possible reasons for this discrepancy: 

(1) the heat budget equations become less accurate as shorter 

periods are considered, 

(2) the change in heat content from one day to another is a 

small difference between two large numbers, 

(3) there may be a loss of cold water through the lateral 

boundaries of the section at depth. 
The third alternative seems to be the most probable. Anderson (195);) 
said that the accuracy of the heat budget for periods greater than 7 days 
Bs + 5 percent. Allowing for the greater accuracy of meteorological data, 
ae heat budget for a period of 22 days should be accurate to at least 
+ 10 percent. The accuracy of the observed heat change also increases 


as longer periods are considered. 
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It is proposed that the computed current be increased by 20 
percent to compensate for the water that flows out of the sides of the 
section at depth. The average flow velocity would then be about 2 cm/sec. 

Assuming an average heat loss of 300 cal/cem?/day, the top 5m of 
the lake would cool at a rate of 0.6° C/day. A flow of 2 cm/sec would 
cause complete replacement of the water in Juanita Bay in 1.5 to 2 days. 
In that time a parcel of water is in the bay would cool 0.9° C to 1.2° C. 
Such cooling corresponds well with the observed temperature decrease in 


Juanita Bay. 
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CONCLUSIONS 


The bottom water in Lake Washington is formed in the many shallow 
bays that surround the lake. The bays temporarily confine the water as 
it cools, allowing it to become denser than the water in the main body 
of the lake. The cold, dense water then flows out of the bay and down 
the slope, mixing with and displacing the original deep water, resulting 
in a bottom layer of nearly uniform temperature. 

The shape of a cross section determines in part whether that section 
has an observed heat loss which is greater or less than the surface heat 
loss computed from the heat budget. Water from sections with extensive 
shallow areas will flow out of those sections into the deep parts of 
adjacent sections that do not have large shallow areas. Continuity 
requires a return flow at some other depth, probably near the surface. 
Because of this removal of cold water from the section in which it was 
formed, sections with large shallow areas will not become as cold as 
expected when heat loss is computed from meteorological conditions. 
Similarly, sections without extensive shallow areas will become colder 
than expected, because of the lateral influx of cold water and loss of 
warmer water. 

Correcting for the escape of cold water through the sides of the 
Juanita Bay section, the estimate of the current velocity was 2 cm/sec. 

Future studies should include cross sections off both shallow 
embayed areas and steep~sided areas off straight shores. They should 
also include routine sections to give a better picture of movements 


along all axes. This however, would require more than one vessel and 
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was not considered within the scope of the present exploratory work. In 


areas where the flow is concentrated, as off the mouth of a bay, a dye 


tracer could be used to follow water parcels and get additional information 


for estimating water velocities. 
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TABLE I 


AVERAGE DAILY WIND FOR THE WARMING CYCLE, 2 FEBRUARY TO 2) APRIL 1963 








Date Wind Date Wind Date Wind Date Wind 
2 Feb. 2.3 23 Feb. 1.5 16 Mar. 1.6 S apr. 2.3 
3 5.6 2h 0.9 i 2.6 6 3.5 
\ 4.3 25 26 18 bs: 7 1.4 
5 Lah 26 5.9 19 Bar 8 on 
6 2.4 27 Lb 20 8 9 lye 
7 ed 28 Diy, 21 Sis 10 2.0 
8 dre 1 Mar. 3.8 22 Bru 11 seg, 
9 2.8 2 2.3 23 2.0 12 251 
10 2.7 3 1.5 2h, as 13 232 
el 1.8 \ Lac 25 Wee 14 3.2 
12 255 5 2.0 26 des 15 6.0 
13 Lats 6 1.9 27 all 16 2.5 
Ly 0.8 i gee 28 ba 17 1.6 
15 @.0 8 1.6 29 Seif 18 oie 
16 0.9 9 1.8 30 5.6 19 2.5 
17 2.2 10 ee an ge 20 1.9 
18 ORs 11 360 1 Apr.) yee? 21 2.4 
19 2.4 12 hes 2 Oey 22 5.1 
20 0.6 13 Gil 3 2.5 23 4.2 
21 1.0 ly 4.2 \ ce. 2h 2a 
22 16S: 15 5.0 
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TABLE II 
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TERMS OF THE HEAT BUDGET EQUATION 


Qr 
22 
22 


20 


Lv) “J Les) —-— Wm 


23 
15 
20 


10 


23 


19 
18 
if 


Fook 


28 


ant 


18 
TABLE III 
Cal/cm*/day 
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FIGURE 12. Top: Daily average air temperature at Sand Point. 
Bottom: Daily average wind speed at Sand Point. 
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Top: Heat loss computed from the heat budget and the observed heat loss in 
the Sand Point - Juanita Bay section. 


Bottom: Computed current in the Sand Point- Juanita Bay sectian. 
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